Observation of spinor dynamics in optically trapped ^'^Rb Bose-Einstein Condensates 
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We measure spin mixing of F=l and F=2 spinor condensates of ^^Rb atoms confined in an 
optical trap. We determine the spin mixing time to be typically less than 600 ms and observe spin 
population oscillations. The equilibrium spin configuration in the F=l manifold is measured for 
different magnetic fields and found to show ferromagnetic behavior for low field gradients. An F=2 
condensate is created by microwave excitation from F=l manifold, and this spin- 2 condensate is 
observed to decay exponentially with time constant 250 ms. Despite the short lifetime in the F=2 
manifold, spin mixing of the condensate is observed within 50 ms. 
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One of the hallmarks of Bose-Einstein condensation 
(BEC) in dilute atomic gases is the relatively weak and 
well-characterized inter-atomic interactions that allow 
quantitative comparison with theory. The vast major- 
ity of experimental work has involved single component 
systems, using magnetic traps confining just one Zee- 
man sub-level in the ground state hyperfine manifold. 
An important frontier in BEC research is the extension 
to multi-component systems, which provides a unique 
opportunity for exploring coupled, interacting quantum 
fluids. In particular, atomic BECs with internal spin 
degrees of freedom offer a new form of coherent matter 
with complex internal quantum structures. The first two- 
component condensate was produced utilizing two hy- 
perfine states in ^^Rb, and remarkable phenomena such 
as phase separation were observed 0, Sodium F=l 
spinor BECs have been created by transferring spin po- 
larized condensates into a far-off resonant optical trap 
to liberate the internal spin degrees of freedom. This 
allowed investigations of the ground state properties of 
Na spinor condensates, and observations of domain struc- 
tures, met ast ability, and quantum spin tunneling 0,0,0. 

In this letter, we explore the spin dynamics and ground 
state properties of ^''Rb spinor condensates in an all- 
optical trap, by starting with well-characterized initial 
conditions in a known magnetic field. We focus on the 
F=l case and confirm the predicted ferromagnetic behav- 
ior. We observe population oscillation between different 
spin states during the spin mixing and observe reduced 
magnetization fluctuations, pointing the way to future 
exploration of the underlying spin squeezing and spin en- 
tanglement predicted for the system 0. We also create 
F=2 spinors using a microwave excitation, measure a de- 
cay of the condensate with a time constant of 250 ms. 
Despite the short lifetime, spin mixing of the spin-2 con- 
densates is observed within 50 ms. Similar results are 
concurrently reported in Ref 7]; in that work, the em- 
phasis is on the F=2 mixing, while here, we focus mainly 
on the F=l manifold. 

A spinor BEC can be described by a multi-component 



order parameter which is invariant under gauge transfor- 
mation and rotation in spin space @, |3, • For a spin- 
1 BEC, the condensate is either ferromagnetic or anti- 
ferromagnetic 0, and the corresponding ground state 
structure and dynamical properties of these two cases 
are very distinct. The Na F=l spinor was found to be 
anti-ferromagnetic, while the F=l ^^Rb is predicted to 
be ferromagnetic [ijj, [l^ • Even richer dynamics are pre- 
dicted for spin-2 condensates [l^, although they remain 
largely unexplored experimentally . 

Single component BECs are typically well-described 
by a scalar order parameter '^(r, t) (the BEC "wave- 
function" ) whose dynamics is governed by the nonlinear 



Gross-Pitaevskii equation p^, ih^ = -^V^^p ^Vt^J ^ 
colV^PV^, where Vt is the trap potential, cq the two-body 
mean field interaction coefficient, and = n is the den- 
sity. For spinor condensates, the Bogoliubov formalism 
can be extended to a vector order parameter ^/^x? where 
X is the spin vector with 2F+1 components for spin-F 
condensates00. For F=l, the two-body interaction en- 
ergy including spin is U{r) = (5(r)(co + C2F1 • F2), where 
r is the distance between two atoms and C2 is the spin 
dependent mean- field interaction coefficient. For F=2, 
U{r) = S{r){a + /?Fi • F2 + 57P0), where a is a spin- 
independent coefficient, /3 and 7 are spin-dependent co- 
efficients, and Pq is the projection operator 00. 

The spin-1 ground states have been calculated by sev- 
eral groups for both zero magnetic field 01 and finite 
field cases 16j. Although the spin dependent terms, e.g. 
|c2 1 , are typically one to two orders of magnitude less than 
the spin independent mean field interaction |co|, they can 
have a dramatic effect on the ground state structure and 
dynamical properties of the condensates. Ho showed that 
if C2 is negative (positive), the spinor displays ferromag- 
netic (anti- ferromagnetic) behavior 0. For the case of 
F=l s^Rb, C2 is calculated to be -3.58(57) x IQ-^^Hz • m^ 
[1^ . Hence the ground state of ^''Rb should be ferromag- 
netic ^3 with a global minimal energy state at zero field 
of |F = l,mF = l) 

For finite fields, the ferromagnetism will manifest it- 



self in different ways depending on the scale of the fer- 
romagnetic energy compared with the Zeeman energy. 
For a typical density of n = 4 x 10^"^ cm~^ in our op- 
tical traps, the ferromagnetic energy 2|c2|n is only 28 
Hz and hence, observation of the low field ground 
state, \F = l,mF = 1), requires that the first order (lin- 
ear) Zeeman shift Ez = tufB x (700 Hz/mG) be smaller 
than 2|c2|n, and hence requires zeroing the magnetic field 
B to less than 40 /iG, typically requiring a magnetically 
shielded environment. 

Even at much larger fields (up to ~500 mG in our case), 
however, the ferromagnetic interactions can still play a 
dominant role in determining the spin ground state due 
to constraints imposed by angular momentum conserva- 
tion. If we start with a non-equilibrium spin mixture, 
the system will relax to the minimal energy state via 
spin exchange collisions, and for the F=l manifold, the 
only processes that conserve angular momentum are 

2|mo)^|mi) + |m_i). (1) 

Hence, the normalized magnetization of the system 
M = ni — n_i is a conserved quantity, where rii is 
the fraction of condensate atoms in the ith Zeeman 
state m^, and for a given M, the system is uniquely 
determined by no- In this field regime, the effect of 
the anti-ferromagnetism (ferromagnetism) is to lower 
(raise) the energy of the mo spinor component relative 
to the average energy of the m±i components 3], which 
drives the reaction in Eq. to the left (right). Hence, 
the extent to which this reaction is driven provides an 
unambiguous distinction between the ferro- and anti- 
ferromagnetic cases. At higher fields, these effects com- 
pete with the second order (quadratic) Zeeman shift, 
{E^i + E_i - 2Eo)/2 B^x (350 Hz/G^) ^^], which 
tends to drive the reaction to the left (to the mo state) 
for fields B > 500 mG. Although evidence of the fer- 
romagnetism of ^^Rb was already provided by the ob- 
served spin mixtures in our previous work (which were 
measured at ^100 mG) [l^^llS, and currently reported 
by Ref. (which are measured at ~340 mG), here we 
present systematic studies of the spin mixing, starting 
with non-equilibrium initial conditions and following its 
time evolution for different magnetic fields. 

Our experiments employ an all-optical BEC technique 
previously described in The atoms are loaded di- 

rectly from a magneto-optical trap into an optical dipole 
force trap formed by a CO2 laser. Lowering the laser 
power forces evaporative cooling in the optical trap, lead- 
ing to rapid condensation. In this work, we employ a 
large period (5 .3 z im) 1-D lattice made by a CO2 laser 
standing wave [l^ which provides a strongly anisotropic 
pancake shape trap, allowing clear distinction between 
thermal clouds and condensates. The lattice is loaded by 
transferring atoms from an orthogonal travelling wave 
trap. We create condensates in only one lattice site by 
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FIG. 1: By applying different magnetic field gradients at dif- 
ferent stages of evaporation, we can create pure condensates 
of a particular spin or mixtures, which consists of 30,000 to 
15,000 atoms as seen above. Each absorptive image is taken 
after 6 ms of free expansion, and a weak Stern- Gerlach field 
is applied during the first 2 ms of expansion to separate three 
spin components spatially. The field of view of each image is 
460 X 230 /im. 
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FIG. 2: Spin mixing of spinor condensate in the optical trap. 
The mixing process starts with a pure ttif = condensate 
in the optical trap. Three separate measurements of the spin 
state are shown after 2 s of spin mixing. 



adjusting the trap power during transfer. The conden- 
sates contain 30, 000 atoms in a single lattice site with 
measured trap frequencies 27r(120, 120, 2550) rad/s, and 
no observable thermal component. The density in the 
optical trap is estimated to be 4.3 x 10^^ cm~^, and the 
Thomas-Fermi condensate radii are (7.6, 7.6, 0.36) /im. 
The 1/e lifetime of the condensates in our optical trap is 
about 3 s. 

To control the initial spin population, we apply differ- 
ent magnetic field gradients during the evaporation pro- 
cess JJJ. To create a pure F = 1, m^ = condensate, we 
apply a field gradient of 28 G/cm during the final 1 s of 
evaporation. To create an equally mixed ttif = —1,0,1 
condensate, a smaller gradient (14 G/cm) is applied. We 
can also create a pure ttif = —1 by applying a 28 G/cm 
gradient at an earlier time before the transfer to the lat- 
tice. Typical spinor condensates with different spin con- 
figurations are illustrated in Figure Q 

To study the spin mixing dynamics, we begin with pure 
rriF = condensates as the initial condition. After con- 
densation, the field gradient is turned off, and a variable 
magnetic bias field is adiabatically ramped up in 10 ms. 
This field can be directed either along the tight (axial) or 
weak (radial) axes of the confinement potential. The con- 
densate is then allowed to evolve for a variable amount 
of time, and then the spin populations are measured us- 
ing absorptive imaging following 6 ms of free expansion. 
To spatially separate the spin components, a weak Stern- 
Gerlach field is applied during the first 2 ms of expansion. 
Typical results for 2 s of spin mixing are shown in Figure 
El We note from these three images taken under identical 
conditions that there is significant variation on the degree 
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FIG. 3: Plot of the fraction in the ttif—O state no and mag- 
netization M (open circle) vs. time. In this measurement, 
M is determined to be 0.5±0.8% initially and -3.5=b2% after 
mixing. Note that there is a 8.5-fold reduction on the fluctu- 
ation of M compare to that of the total population (which is 
measured to be 17%). Clear population oscillations of no are 
seen, and the fluctuation in no is 6.6 fold of that of M. This 
data was taken at a bias field of 100 mG 

of mixing from run to run of the experiment. However, 
in each case, the magnetization appears to be conserved. 
Generally, the components of the spin mixtures are iden- 
tical in shape to the original cloud to within our imaging 
resolution. Occasionally one or more of the components 
will appear to have either a thermal component or a dis- 
torted shape. 

We have measured the spin mixing for different evolu- 
tion times following preparation of the rriF = conden- 
sates. The time to reach equilibrium of the spin mixing 
is typically less than 600 ms, and this time decreases 
slightly with increasing magnetic field. Fig. [2l shows the 
average time evolution of no and M for 50 repeated mea- 
surements. First, we note that the magnetization, M, is 
conserved throughout the mixing to within a few percent. 
Although the data does suggest a drift of M below zero 
by a small amount -3.5±2% (the uncertainty here and er- 
ror bars in the data are purely statistical), the deviation, 
if any, is comparable to our uncertainties in measuring 
populations (~3%), limited by the absorptive imaging 
technique. Secondly, there is an almost 10-fold reduction 
in the statistical noise of the magnetization relative to 
that of the total population, which varies 17% from shot 
to shot. This suggests that the fluctuations of ni and 
n_i, which are coupled from the mixing processes in Eq. 
(QJ, are quantum correlated. These correlations underly 
theoretical predictions for spin squeezing and entangle- 
ment in the system j^ij. Thirdly, the relaxation of no 
to the steady state value is not monotonic but instead 
shows a damped oscillation at 4 Hz. Such oscillations 
are a natural outcome of coherent spinor mixing as shown 
theoretically in glQj . 
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FIG. 4: Plot of no vs. magnetic field after 3 s of spin mix- 
ing. The upper curves are the theoretical predictions for the 
ant i- ferromagnetic case, while the lower curves are for the fer- 
romagnetic case. The dashed lines are the predictions with 
a field of gradient 20 mG/cm applied, which is the measured 
upper bound in our trap. 



We attempted to directly measure the phase relation- 
ship between the spinor components by performing 
an interference experiment between two spinor conden- 
sates created in adjacent lattice sites. In the absence of 
dephasing mechanisms, we would expect to see clear in- 
terference fringes in time-of-flight measurements both be- 
fore and after spin mixing. In the measurement, clear in- 
terference fringes were visible initially, while fringes were 
typically not observed after mixing. It is quite possible 
however that the observed dephasing of the interference 
pattern was caused by small field variations between two 
sites due to a stray magnetic field gradient (< 20 mG/cm) 

To make comparison with theoretical predictions 
we measure the degree of mixing for different applied 
magnetic fields. Fig. ^ shows the results of spin mix- 
ing for 3 s, in which no is plotted vs. the applied field. 
Also shown in this figure are theoretical curves for both 
the antiferro- and ferromagnetic cases in a magnetic field 
with and without a small (20 mG/cm) field gradient 16]. 
As evidenced by the data, the spin mixing agrees with the 
ferromagnetic predictions and is inconsistent with the 
anti-ferromagnetic prediction. When the field is larger 
than 700 mG, the quadratic Zeeman effect completely 
dominates the spin interaction, and the condensates re- 
main in the = state. Note that magnetic fields are 
applied either along tight trap or weak trap axes with 
different polarities; however, no significant difference in 
the measurements is observed. 

We also measured spin mixing for fields less than 100 
mG but found that our results were affected by the 
stray ~10 mG AC magnetic fields present in the cham- 
ber. These fields are capable of driving off-resonant rf- 
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FIG. 5: (a) Lifetime measurement of F=2 spinor condensates. 
Following a rapid decay initial decay, the population decays 
exponentially with a time constant of 250 ms. (b) Spin mixing 
of the F=2 spinor condensates in the optical trap. These 
images represent three identical measurements after 50 ms of 
spin mixing. Note that the magnetization is conserved. 



transitions between the Zeeman sub-levels. We observed 
this directly by creating aF = l,mF = — 1 condensate 
(for which there is no spin mixing due to conservation 
of magnetization) and measuring the final spin popu- 
lation. For magnetic fields greater than 100 mG, the 
magnetization remains conserved, while at lower fields, 
the niF = —1 atoms are quickly pumped by the AG 
fields (within 100 ms for fields <10 mG) to other Zee- 
man states. 

To study F=2 spinors, we coherently excite the pure 
F = l,mi? = condensates to F = 2,mi? = using 
microwave fields tuned to 6.8 GHz. Additionally, by con- 
trolling the bias field, the microwave frequency, and the 
initial Zeeman state in the F=l manifold, we can pump 
the condensate to any Zeeman sub-level of the excited hy- 
perfine manifold. The F=2, m=0 condensate is observed 
to decay as shown in Fig [Sfa). Following an initially 
rapid decay, it decays exponentially with a time constant 
of 250 ms. Despite the short lifetime in the excited hyper- 
fine manifold, we still observe spin mixing within 50 ms, 
and magnetization conservation is also observed during 
the mixing, as shown in Fig.jSfb). 

In summary, we have observed spin mixing of ^^Rb 
spinor condensates in F=l and F=2 hyperfine mani- 
folds in an optical trap. The observed equilibrium spinor 
configurations of the lower manifold confirms that F=l 
^^Rb is ferromagnetic. The magnetization was conserved 



within the measurement errors during the entire spin 
mixing process. The reduced noise of magnetization sug- 
gests quantum correlation of the spin dynamics, which 
underlies spin squeezing and spin entanglement. Future 
work will study the coherence of spin mixing, spin squeez- 
ing, and entanglement. 
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